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Abstract 

Bose-Einstein correlations of pairs of identical charged pions produced 
in hadronic Z decays are analyzed in terms of various parametrizations. 
' A good description is achieved using a Levy stable distribution in con- 

junction with a model where a particle's momentum is highly correlated 
with its space-time point of production, the r-model. However, a small 
' but significant elongation of the particle emission region is observed in 

the Longitudinal Center of Mass frame, which is not accommodated in 
the r-model. This is investigated using an ad hoc modification of the 
0^ ' r-model. 
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We study Bose-Einstein correlations (BEC) in hadronic Z decay using data 
(— I . collected by the l3 detector at an e'^'e" center-of-mass energy of ^/s ~ 91.2 

GeV. Approximately 36 million like-sign pairs of well-measured charged tracks 
from about 0.8 million hadronic Z decays are used [1]. Events are classified as 
two- or three-jet events using calorimeter clusters with the Durham jet algorithm 
. with a jet resolution parameter j/cut = 0.006, yielding about 0.5 million two-jet 

\^ ' events. There are few events with more than three jets, and they are included 

, in the three-jet sample. To determine the thrust axis of the event we also use 

calorimeter clusters. 

lO . The two-particle correlation function of two particles with four-momenta pi 

' and p2 is given by the ratio of the two-particle number density, p2(pi,P2), to the 

product of the two single-particle densities, pi (pi)pi {pt)- Since we are interested 
only in BEC, Pi{pi)pi{p2) is replaced by pq{pi,P2), the two-particle density that 
would occur in the absence of BEC, yielding i?2(pi,P2) = P2{pi,P2)/ Pa{pi,P2)- 
' An event mixing technique is used to construct po- 

^ , Since the mass of the identical particles of the pair is fixed, i?2 is defined in 

six-dimensional momentum space, which is often reduced to a single dimension, 
Q — \/—{pi — P2)^- But there is no reason to expect the hadron source to be 
spherically symmetric in jet fragmentation. In fact, the source is found to be 
elongated along the jet axis [^-[5], but only by about 25%, which suggests that a 
parametrization in terms of the single variable Q, may be a good approximation. 
This is also suggested by studies of various decompositions of Q [Tl|6] . 
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2 Parametrizations of BEC 



With a few assumptions, R2 is related to the Fourier transform, f{Q), of the 
(configuration space) density distribution of the source, f{x): 



R2{PuP2) = l l + A|/(Q)p {1 + 5Q) 



(1) 



The parameter 7 and the (1 + 5Q) term are introduced to parametrize possible 
long-range correlations not adequately accounted for in poj £ind A to account for 
several factors, such as lack of complete incoherence of particle production and 
presence of long-lived resonance decays. 

Model-independent ways to study deviations from the Gaussian are to ex- 
pand f{Q) about a Gaussian using the Edgeworth expansion or to replace f{Q) 
by a symmetric Levy distribution. 

Both the symmetric Levy (Fig. la) and the Edgeworth parametrizations do 
a fair job [Zl[8] of describing the region Q < 0.6 GeV, but fail at higher Q, 
particularly the region 0.6-1.5 GeV where R2 dips below unity, indicative of an 
anti-correlation. This is clearly seen in Fig. la. 
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Figure 1: The Bose- Einstein correlation function R2 for two-jet events. The 
curves correspond to fits of the parametrizations: (a) the symmetric Lcvy;(b) 
the simplified r- model, Eq. (|4]); and (c) the full r- model, Eq. ([2]). Also plotted 
is A, the difference between the fit and the data. The dashed line represents 
the long-range part of the fit, i.e., 7(1 + eQ). The dotted line in (a) represents 
a linear fit in the region Q > 1.5 GeV. 

We have seen that BEC depend, at least approximately, only on Q, not on 
its components separately and that R2 in the region 0.6-1.5 GeV dips below 
its values at higher Q. A model which predicts such Q-dependence as well 
as a dependence on the transverse mass, mt = \/ m? + , which has been 
previously observed [5HTT| , is the r- model [12] , in which it is assumed that in the 
overall center-of-mass system the average production point x = (t,fx,fy,rz), of 
particles with a given four- momentum p is given by x^{jp^) = arp^ . In the case 

of two-jet events, a = 1/mt, and r = \Ji^ — f\ is the longitudinal proper time; 
for the case of three-jet events the relation is more complicated. The second 
assumption is that the distribution of x^(j)^^ about its average is narrower than 
the proper-time distribution, H{t). Then R2 is found |13] to depend only on 
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Q, the values of a of the two pions, and the Fourier transform of H{t). Since 
there is no particle production before the onset of the collision, H{t) should 
be a one-sided distribution. We choose a one-sided Levy distribution, which is 
characterized by three parameters: the index of stability a, the proper time of 
the start of particle emission tq, and At, which is a measure of the width of 
H{t). Then [B] 



i?2(Q, 01,02) = 7 < H- Acos 



ToQ^(ai -I- 02) 



tan 



• exp 



(1 + eQ) . 



(2) 



Note that the cosine factor generates oscillations corresponding to alternat- 
ing correlated and anti-correlated regions, a feature clearly seen in the data 
(Fig. 1). Note also that since a = l/irit for two-jet events, the r-modcl predicts 
a decreasing cfFcctivc source size with increasing mt. 

Before proceeding to fits of Eq. we first consider a simplification of the 
equation obtained by assuming (a) that particle production starts immediately, 
i.e.. To — 0; and (b) an average a-depcndence, which is implemented by intro- 
ducing an effective radius, i?, defined by 



(3) 



This results in 



i?2(Q)=7 1 + Acos ((i?aO)'") exp (-(i?gf") (1 + eQ), (4) 



where i?a is related to R by 



R. 



2a 



.a.,(^)«- 



(5) 



Fits of Eq. Q are first performed with i?a as a free parameter. The fits for both 
two- and three-jet events have acceptable confidence levels (CL), and describe 
well the dip in the 0.6-1.5 GeV region, as well as the peak at low values of Q. 
The estimates of some fit parameters arc rather highly correlated. For example, 
for two-jet events the estimated correlation coefficients from the fit for a, R 
and i?a are p{a,R) = -0.62, p{a,R^) ^ -0.92, and p{R,R^) ^ 0.38. Taking 
the correlations into account, the fit parameters satisfy Eq. ([5]), the difference 
between the left- and right-hand sides of the equation being less than 1 standard 
deviation for two-jet events and about 1.5 standard deviations for three-jet 
events. No significant long-range correlation is observed: e is zero within 1 
standard deviation, and fits with e fixed to zero find the same values, within 1 
standard deviation, of the parameters as the fits with e free. Thus the method 
to remove non-Bose-Einstein correlations in the reference sample is apparently 
successful. 

Fits are also performed imposing Eq. ([S]). The result for two-jet events is 
shown in Fig. lb. For two-jet events, the values of the parameters arc compara- 
ble to those with i?a free. For three-jet events, the imposition of Eq. ([5]) results 
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in values of a and R closer to those for two-jet events, but the is notice- 
ably worse, though acceptable, than with i?a free. Also, the values of e differ 
somewhat from zero, which could indicate a slight deficiency in the description 
of BEC similar to, but on a much smaller scale, the Edgeworth and symmetric 
Levy fits. 

For two-jet events, a = l/wt, while for three-jet events the situation is more 
complicated. We therefore limit fits of Eq. ^ to the two-jet data. For each bin 
in Q the average values of mti and mt2 are calculated, where mti and TOt2 are the 
transverse masses of the two particles making up a pair, requiring mti > mt2- 
Using these averages, Eq. (0) is fit to i?2(Q)- The fit yields tq = 0.00 ± 0.02 fm, 
and we have re-fit with tq fixed to zero. The results are shown in Fig. Ic. 
The parameters a, Ar and A are highly correlated with estimated correlaton 
coeflicients p[a, At) = -0.96, p[a, A) = -0.92, and /9(At, A) = 0.95. 

Since the r-model describes the rat dependence of R2, its parameters, a, Ar, 
and To, should not depend on Wf However, A, which is not a parameter of the 
r-model, but rather a measure of the strength of the BEC, can depend on mt- 
The large correlation between the fit estimates of A, a, and Ar complicate the 
testing of mt-independence. We perform fits in various regions of the mti-mt2 
plane keeping a and Ar fixed at the values obtained in the fit to the entire 
TOt plane. The regions are chosen such that the numbers of pairs of particles 
in the regions are comparable. The mt regions and the CLs of the fits are 
shown in Fig. (5] The CLs are reasonably uniformly distributed between and 
1. The data are thus in agreement with the hypothesis of mt-independence of 
the parameters of the r-model. 

3 Do BEC depend on components of Q? 

In this section we return to the question of 
a possible elongation of the region of homo- 
geneity, as viewed in the rest frame of the pion 
pair, and the possible dependence of the BEC 
correlation function on components of Q. 

The r-model predicts that the two- 
particle BEC correlation function R2 depends 
on the two-particle momentum difference only 
through Q, not through components of Q. 
However, previous studies [2 H51ITT] have found 
that i?2 does depend on components of Q 
and indeed that the shape of the region of 
homogeneity is elongated along the event 
(thrust) axis. The question is whether this 
is an artifact of the Edgeworth or Caus- 
sian parametrizations used in these studies or 
shows a defect of the r-modcl. 

In the studies of elongation, is split 
into three components in the Longitudinal 
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Figure 2: The regions of the 
mti-mt2 plane and the CLs of 
the fits in these regions. 
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Center of Mass System (LCMS) of the pair: 



Q=Qi + Osidc + Qout - i^EY (6) 

= Ql + QLe + QLt (1 - P') , where /? = ^^g+g^" . (7) 

The LCMS frame is defined as the frame, obtained by a Lorcntz boost along the 
event axis, where the sum of the three- momenta of the two pions (pi + P2) is 
perpendicular to the event axis. Assuming azimuthal symmetry about the event 
axis suggests that the region of homogeneity have an ellipsoidal shape with one 
axis, the longitudinal axis, along the event axis. 

In the LCMS frame the event axis is referred to as the longitudinal direction, 
the direction of pi + p2 as the out direction, and the direction perpendicular to 
these as the side direction. The components of \pi — p2 \ along these directions 
are denoted Ql, Qout, and Qsidc, respectively. R^Q^ is then replaced in the 
Gaussian or Edgeworth parametrizations by 

i?2Q2 ^ RIQI + i?sidoQsidc + Pout Qout ■ (8) 

i?L and i?sidc are measures of the longitudinal and transverse size of the source, 
respectively, whereas the value of pout reflects both the transverse and temporal 
sizes H 

Since the present analysis uses a somewhat different event sample and mixing 
algorithm, we have first checked that we find values of Rsidc/ Ri^ consistent with 
our previous analysis [2. 

We next investigate whether the r-modcl parametrizations could accommo- 
date an elongation by incorporating Eq. (|5]) in Eq. If Eq. ([5|) is imposed, this 
results in 



1 + A cos (tan A^" ) exp {~A^ 



(1 + eLQL + EsidcQsidc + EoutQcut) 

(9) 



i?2(Q)=7 
where 

= RlQl + i?LoQLe + Pout Qout • (10) 

A fit of Eq. ^ for two-jet events for Q < 4 GeV has a CL of 66%, and a clear 
preference for elongation is found: i?side/^L = 0.61 ± 0.02. The value of R 
found in the r-model fit shown in Fig. Ic lies between the values of i?sido and 
i?L found here. Furthermore, the value of Rgide/Rh agrees well with that found 
using the Edgeworth parametrization. 

To directly test the hypothesis of no separate dependence on components of 
Q, we investigate two decompositions of Q: 

Q' = Qle + QLc + Qout , Qle = Ql - i^E f ; (11a) 

Q' = QL + QsMc + '?out, <z,\t = - (Ai?)^ (lib) 

The first, Eq. (jllap . corresponds to the LCMS frame where the longitudinal and 
energy terms are combined; its three components of Q are invariant with respect 



to Lorentz boosts along the event axis. The second, Eq. (jllb[) . corresponds to 



^In the literature the coefficient of Q^^j in Eq. ^ is usually denoted R^^^^- We prefer to 
use to emphasize that, unlike -Rl and -Rsido Pout contains a dependence on ^. 
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the LCMS frame boosted to the rest frame of the pair; its three components are 
invariant under Lorentz boosts along the out direction. 

The test then consists of replacing E?Q^ by times one of the above 
equations and comparing a fit where the coefficients of all three terms are con- 
strained to be equal to a fit where each coefficient is a free parameter. For both 
parametrizations the fit allowing separate dependence on the components of Q 
has a significantly better CL than the fit which does not. The fit using Eq. (jllap 
attains a CL of 2% when no separate dependence is allowed. While this is not 
poor enough to reject, by itself, the hypothesis of separate dependence, the fit 
allowing it is, with CL = 38%, much better. When Eq. (jlfbp is used the fit 
not allowing separate dependence is rejected by a CL of 10"*", while allowing 
it results in the marginally acceptable CL of 2%. For both decompositions the 
differences in between the fit allowing separate dependence and that not 
allowing it is huge, 296 and 464, respectively, while a difference of only 2 is 
expected if there is no separate dependence. This provides very strong evidence 
in favor of separate dependence on the components of Q. Moreover, the radii 
corresponding to Qsidc and gout are found to be unequal by about a factor 2, 
which violates the assumption of azimuthal symmetry assumed in the LCMS 
analysis. 

4 Discussion 

The usual parametrizations of BEC of pion pairs in terms of Q, such as the 
Gaussian, Edgcworth, or symmetric Levy parametrizations, are found to be 
incapable of describing the BEC correlation function i?2, particularly in the 
anti-correlation region, 0.5 < Q < 1-5 GeV. 

The T-model assumes a very high degree of correlation between momentum 
space and space-time and introduces a time dependence explicitly. The two- 
particle BEC correlation function i?2 is then dependent on one two-particle 
variable, Q, and on a, a parameter in the correlation between momentum space 
and space-time. The time-dependence is assumed to be given by an asymmetric 
Levy distribution, H{t). The parameter a can be absorbed into an effective 
radius to obtain Eq. which successfully describes BEC for both two- and 
three-jet events. 

For two-jet events a = 1/mt and the introduction of an effective radius is 
unnecessary. The BEC correlation function is then a function not only of Q but 
also of the transverse masses of the pions. This description of R2{Q, rriti, mt2) 
also successfully describes the data. 

Nevertheless, the r-model description breaks down when confronted with 
data expressed in components of Q. The r-model predicts that the only two- 
particle variable entering R2 is Q. There is no dependence on components of 
Q separately. This implies a spherical shape of the pion region of homogene- 
ity, whereas previous analyses using static parametrizations have found a shape 
elongated along the event axis, as evidenced by i?side being smaller than i?L in 
the LCMS. Assuming azimuthal symmetry about the event axis, i?sido is the 
transverse radius of the region of homogeneity. Accordingly, the r-modcl equa- 
tions for i?2 have been modified ad hoc to allow an elongation. The elongation 
found agrees with that found in previous analyses. However, in the rest frame of 
the pair the radius parameter in the out direction is found to be approximately 
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twice that in the side direction. Thus the assumption of azimuthal symmetry 
about the event axis is found not to hold. 

We note that elongation is expected [l3] for a hadronizing string in the 
Lund model, where the transverse and longitudinal components of a particle's 
momentum arise from different mechanisms. Perhaps the absence of azimuthal 
symmetry arises from gluon radiation, which is represented by kinks in the Lund 
string. Therefore a possible improvement of the r-model could be to allow the 
longitudinal and transverse components to have different degrees of correlation 
between a particle's momentum and its space-time production point. 
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